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We present results of an improved analysis of the experimental data of the EAS Cherenkov array Tunka-25. 
A new function to fit the Cherenkov light lateral distribution LDF at core distances from to 350 m has been 
developed on the base of CORSIKA simulations and applied to the analysis of Tunka data. Two methods to 
estimate the EAS maximum position have been used. The one is based on the pulse FWHM, the other on the 
light LDF. We present the primary energy spectrum in the energy range 10 15 - 10 17 eV. The use of the depth of 
the EAS maximum to determine the mean mass composition is discussed. 



Tunka-25: Experiment and Simulatons 

The Tunka EAS Cherenkov array is located in Tunka Valley, at an altitude of 675 m a.s.l., and was described 
inRef. Q. 

Extensive air showers in an energy range 3 ■ 10 14 to 2 ■ 10 16 eV have been simulated with CORSIKA. The total 
amount of 600 CORSIKA events has been simulated for primary protons and iron nuclei and for three zenith 
angles 0°, 15° and 25°. The LDF consists of two branches: one almost exponential from the core to some 
distance R kn , the other following a power law from Rk n up to 350 m: 

n( m _{ Q kn -exp{{R kn -R)-(l + 3/{R + 3))/Ro), forR<R kn 

^ ' 1 Qkn ■ {Rkn/R)\ forR>R kn W 

i?o = 10 2 - 95 - a245P ,m 
Rkn = 155 - 13P, m 
b= 1.19 + 0. 23P 

The main difference to the expression suggested in 1 2 1 is the use of variable power law index for the second 
branch. All three partial parameters defining the LDF shape, Rq, R kn and b, are strict functions of a single 
steepness parameter P, defined as the ratio of Cherekov light fluxes at core distances 100 and 200 m: P = 

Q(100)/Q(200). 

Examples of fitting the simulated LDF with equation 1 for different parameters P are shown in fig. 1 . The 
parameter P and the light flux at a fixed core distance 175 m, Q175, have been determined for every simulated 
event. 

The essential correlations of parameter pairs E n and Q175, Eq and H max , H max and P, and the standard 
deviations of their distibutions, separately for p and Fe, have been extracted from this CORSIKA simulation 
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Figure 1. CORSIKA: EAS Cherenkov light LDFs and fit- Figure 2. EXPERIMENT: Differencial energy spectrum, 
ting functions. 1 - P = 5, 2 - P = 4, 3 - P = 3 



in order to use them in a special code called "model of experiment". Two intermediate nuclear groups, cor- 
responding to primaries He and CNO have been added to the "model" using an interpolation of parameters, 
obtained for p and Fe, The real geometry of the array and the response of every detector and it's fluctuations 
are taken into account in the model. The "model of experiment" permits a fast generation of 10 4 — 10 5 ar- 
tificial events, distributed with a realistic energy spectrum, for different assumptions concerning the primary 
mass composition. It allows to analyse the array response, including the analysis of details of the program of 
EAS parameter reconstruction, the efficiency of registration of EAS, the threshold of the total array, errors and 
possible distortions of the reconstructed energy Eq, and the depth of the shower maximum X max . 

The "model of experiment", assuming a complex mass composition (p:He:CNO:Fe=0. 25:0. 25:0. 25:0. 25), has 
shown that for energies > 3 ■ 10 15 eV, our procedures of relative calibration of the detectors and the EAS 
parameter reconstruction result in accuracies of primary energy determination of ~ 15%, of the core position 
of ~5 m and of X max of ^30 g/cm 2 . For energies below 3 • 10 15 the errors are larger. The threshold of data 
acquisition with 100% (50%) efficiency is 8 • 10 14 eV (5 • 10 14 eV). Some other conclusions of the "model" 
will be discussed later. 

Primary Energy Spectrum 

The primary energy Eq [TeV] has been obtained from Q175 [photon • cm~ 2 • eV~ 4 ] with the relation: 

E =400-Q° 17 f. 

The absolute energy calibration is based on the results obtained with the QUEST experiment 1 3 1. 

Figure 2 presents the differential energy spectrum, derived from data taken in 300 hours, spread over about 50 
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clear, moonless nights, with a trigger rate above 1.8 Hz . To construct a spectrum, showers with zenith angles 
9 < 25° and a core position inside the geometrical area of the array have been selected. 

We note that particularly at the lowest and the highest energies, the new fit functions performs significantly 
better than that used in |4 1. 



Depth of the EAS Maximum 

The lateral and the time distributions of the Cherenkov light provide two independent methods to estimate the 
depth of the EAS maximum. The first is the measurement of the LDF steepness P, which is related to the 
distance to the shower maximum by the expression H max (in [km]): H max — 17.63 — 0.0786 * (P + 8.916) 2 . 
This relation is almost independent of other details of the simulation: energy, sort of nucleus, zenith angle and 
model of hadron interaction. The distribution of H max for a fixed P has a standard deviation of only 0.3 km. 

The depth of the EAS maximum X max is derived from H max using the atmosphere parametrs consistent with 
really observed (the mean temperature during the observations is -20°C). The "model of experiment" gives an 
error of 30 g/cm 2 for the experimental depth of the EAS maximum X max . 

The second method is measuring the Cherenkov pulse full width on half maximum (FWHM). In accordance 
with CORSIKA simulations as described above, the FWHM [ns] at distances larger than 200 m from the EAS 
axis is related to the relative position of the EAS maximum: AX = Xq/ cos 9 — X max [g/cm 2 ], - where Xo 
is the total depth of the atmosphere and 9 is the zenith angle of the shower. The relation between FWHM 
and AX depends only on the distance to the EAS axis, and is almost independent on the other details of the 
simulation: energy, sort of nucleus and model of hadron interaction. For example, for a distance of 250 m, one 
obtains AX = 1677 + 1006 • log 10 (FWHM). The underlying theoretical uncertainties for this method are 
smaller than for the first one. 

Figure 3 presents the mean depth of the EAS maximum, derived with the two methods described above, as a 
function of primary energy. The "model of experiment" shows that there is no influence of the array threshold 
on the obtained X max , starting from an energy 1.5 • 10 15 eV. It is seen from fig. 3 that the threshold of the 
FWHM method is higher than that of the LDF steepness method, but the mean depths, obtained with the two 
different methods are in good agreement. Both methods result in large fluctuations of mean points at low 
statistics (high energy), as expected for a very asymmetric distribution like that of X max . 

Experimental X niax distributions for EAS in the energy range 3-10 PeV are shown in fig 4. All the values of 
X max are re-normalized to a fixed energy of 5-10 15 eV, using an experimental elongation rate of (97±3) g/cm 2 . 
The experimental distribution is compared to distributions simulated with the "model of experiment", for dif- 
ferent assumptions on primary composition, and re-normalized in the same way as the experimental one. The 
agreement with experimental data is better for the complex composition (p:He:CNO:Fe=0. 25:0. 25:0. 25:0. 25), 
than for pure p or Fe. So the intermediate mass composition is more probable, but it should be noticed that the 
difference between simulated and experimental standard deviation is significant only for pure p or Fe. For a 
wide range of intermediate compositions, mixtures of mass components lead to standard deviations similar to 
those marked as "complex" in the figure. 

The mean value of the simulated distribution depends on hadron interaction model. In fig. 3, -X" max vs. energy 
is displayed for two models QGSJET-01 and QGSJET-II 0. All other models used in CORSIKA yield 
lines between these borders |[7). 

One sees, that a conclusion about the mean mass composition strongly depends on the chosen hadron interac- 
tion model. To solve the problem with the model the measurement of energy dependence of X max in much 
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Figure 3. Mean depth of EAS maximum vs energy. 



Figure 4. X max distribution. 



wider energy range is very essential. This is one of the reasons of our decision to spread our measurement to 
lower energy with the new version of Tunka-25 Cherenkov light detectors supplied with spheric mirrors of 1 
m diameter. To the other hand the new experiment Tunka-133 1 8 1 with effective area 10 times higher than that 
of Tunka-25 is under construction in Tunka valley. 
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